Controlling the second-harmonic in a phase matched negative- index metamaterial 
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Nonlinear metamaterials (NLMMs) have been predicted to support new and exciting domains in 
the manipulation of light, including novel phase matching schemes for wave mixing. Most notable is 
the so-called nonlinear-optical mirror, in which a nonlinear negative- index medium emits the gener- 
ated frequency towards the source of the pump. For the first time, we experimentally demonstrate 
the nonlinear-optical mirror effect in a bulk negative-index NLMM, along with two other novel 
phase matching configurations, utilizing periodic poling to switch between the three phase matching 
domains. 



Nonlinear metamaterials (NLMMs) represent the pin- 
nacle of designer electromagnetic materials, allowing one 
to manipulate light based on direction, polarization, fre- 
quency, and intensity. In particular, these artificial medi- 
ums have the potential to make nonlinear-optical phe- 
nomena accessible at low powers and in novel and excit- 
ing configurations [TJ [2] . While a wide range of nonlinear 
processes have already been demonstrated in NLMMs, 
including frequency generation [3HS] , parametric amplifi- 
cation [6j [7] , and bistability [8] [9] , the setups have been 
constrained to sub- wavelength interaction lengths for fear 
of the destructive effects of phase mismatch. As such, 
the conversion efficiencies achieved are only a fraction of 
their potential. However, phase matching is an exciting 
subject that involves the entire set of linear and non- 
linear properties, bringing to bear the full capabilities 
of NLMMs to manipulate light at will. This claim has 
been supported by several theoretical studies on NLMMs, 
most notably the prediction of a nonlinear-optical mirror 
that generates second-harmonic (SH) back towards the 
fundamental frequency (FF) source [6j[10j[TT], but exper- 
imental verification of these novel phase matching effects 
in a bulk NLMM has been lacking. 

In this letter, we experimentally demonstrate three 
unique phase matching configurations for second- 
harmonic generation (SHG) in a metallic waveguide 
loaded with a magnetic NLMM at microwave frequencies. 
The three configurations are reflected SH phase match- 
ing in a negative-index band, transmitted SH quasiphase 
matching (QPM), and simultaneous QPM of both the re- 
flected and transmitted SHs near a zero-index band. The 
resulting SH spectrums, showing strong phase matching- 
induced enhancements, are supported qualitatively by 
nonlinear transfer matrix method (TMM) calculations on 
a homogeneous slab with equivalent effective properties. 

The NLMM under investigation in this letter is 
the varactor-loaded split-ring resonator (VLSRR), the 
canonical subject of numerous recent studies. The con- 
stituent unit-cell of this NLMM is composed of a 17 /xm 
thick, 0.75 mm wide copper ring printed on a 0.2 mm FR4 
PCB substrate. The 1 mm gap is loaded with a Sky works 
SMV1231 varactor diode, whose capacitance varies with 
applied voltage according to C(Vd) — Co(l — Vd/Vp) M \ 



where Vd is the voltage across the diode, Co = 2.35 pF is 
the zero- volt age capacitance, V p = 1.5 V is the intrinsic 
potential, and M = 0.8 is the gradient coefficient [T2] . 
Our bulk material, displayed in Fig. [I] consists of 16 lay- 
ers (24 cm in the direction of propagation) of VLSRRs 
placed in an aluminum waveguide with cross-section 15 
cm x 1.5 cm, completely filling the waveguide in the 
transverse direction, and arranged in a square lattice with 
1.5 cm spacings. For the purpose of periodic poling, the 
VLSRR sample is divided into four identical sections that 
can be reoriented independently. The particular choice 
of dimensions ensures that the resonance frequency of 
the VLSRR is below the cutoff frequency of the waveg- 
uide's lowest-order mode. As such, the NLMM loaded 
waveguide is expected to support a backward mode at 
frequencies where the VLSRR's permeability is negative, 
that is, a mode whose phase velocity is opposite to the di- 
rection of power flow [T3HT5] . In addition, near the cutoff 
frequency, the wave vector of this mode vanishes, yielding 
a second region of interest for phase matching [16] . 




FIG. 1: (Color online) Photograph of the waveguide with the 
top removed, loaded with four identical sections of VLSRRs. 
The inset shows an enlarged view of the NLMM unit-cell. 

The linear properties of the VLSRR medium were char- 
acterized in a transmission line setup by measuring the 
transmission of a single layer and fitting to known prop- 
erty models [17], shown in Fig. [2] (a). We find our VL- 
SRR medium to be well-described by a constant permit- 
tivity c y (u)) = 2.2eo, and a permeability given by the 
Lorentz oscillator-like formula 
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where F = 0.22 is the oscillator strength, cuq = 2tt x 608 
MHz is the resonance frequency, and 7 = 2ir x 14 MHz 
is the damping coefficient. The permeability is assumed 
to be that of free-space along the other axes. 
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FIG. 2: (Color online) Plot of the magnitude of the retrieved 
(solid blue) and fitted (dashed green) ^-parameters (a) and 
second-order susceptibility (b) for the VLSRR medium. 

Subsequently, we measured the SHG output from this 
setup. Using the nonlinear TMM retrieval method [T7] - 
[19], we extracted the second-order magnetic nonlinear 
susceptibility from the SH spectrum (shown in Fig. [2] 
(b)). The result was then fitted to an effective medium 
model obtained via a perturbative circuit analysis, [20 
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where a = ±0.198 V -1 is the fitted strength of the non- 
linearity, A = 113 mm 2 is the ring's area, and D(uo) = 
(uJq — ijuj — uj 2 ) is the Lorentz denominator. The sign of 
a is determined by the orientation of the varactor diode 
in the VLSRR, reflecting the fact that reversing the diode 
direction is analogous to reversing the orientation of the 
nonlinear metacrystal, and thus is expected, by symme- 
try arguments, to induce a tt shift in the phase of the 
even-order nonlinear it ies. Note that the magnitude of a 
is close to the predicted value of ^- = 0.267 V" 1 [20]. 

The dispersion of the TE10 mode of the metallic waveg- 
uide is modeled by 



k z (uj) = ±0J\Ui x 
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where k z (uo) is the longitudinal wavevector, uo c — 
j—^= ^ = 2tt x 674 MHz is the cutoff frequency, and 
b = 15 cm is the waveguide width. As in typical negative- 
index media, the negative sign is chosen when both terms 
in the square root are negative, i.e., when ji x < and 
€ y — ^ < 0. The former occurs just above the 

resonance frequency of the VLSRR, while the latter re- 
quires operation of the waveguide below cutoff, so that 
we can expect a backward- wave to exist for frequencies of 
roughly 615 MHz < uj/2tt < 674 MHz. Furthermore, for 



a FF below cutoff, it can be shown that only the TE10 
mode will propagate at the SH frequency, simplifying our 
analysis. For coupling between the TE10 modes at the 
FF and SH, we obtain an overlap integral of T « 1.2 [21]. 

The phase matching effects in our VLSRR loaded 
waveguide can be examined by employing the above ex- 
perimentally fitted models to calculate the SHG phase 
mismatch and coherence lengths. For example, the phase 
mismatch for SHG in a homogeneous medium is given by 



Ak = ±k z (2u) - 2k z (u), 



(4) 



where it is assumed that the energy of the FF wave prop- 
agates in the positive z direction, and the ± refers to a 
positively (transmitted) or negatively (reflected) prop- 
agating SH wave. Thus, assuming the SH is far from 
resonance, it is clear that a forward- wave (k z (u) > 0) 
at the FF will tend to phase match with a transmit- 
ted SH wave, as is seen in conventional nonlinear optics. 
A backward- wave (k z (uj) < 0) at the FF, on the other 
hand, will tend to phase match with a reflected SH wave, 
as predicted earlier [6j [10] [11] . This odd behavior can be 
seen more clearly by computing the coherence lengths, 
L coh = shown in Fig. [3] using Eqs. ^ and Q. As 
expected, the coherence length for the transmitted SH is 
small (less than the interaction length), while the coher- 
ence length for the reflected SH is much larger throughout 
the negative-index band. The reflected phase-matched 
region, shown in gray in Fig. [3] represents the frequency 
range in which our NLMM acts as the nonlinear-optical 
mirror of Refs. [6 and [11 , with the vast majority of the 
SH power traveling opposite to the FF. 









transmitted SHG 








reflected SHG 















1100 1150 1200 1250 1300 1350 1400 

SH frequency (MHz) 

FIG. 3: (Color online) Plot of the calculated coherence lengths 
for transmitted (blue) and reflected (green) SHG in the VL- 
SRR loaded waveguide. The dashed lines at 12 and 24 cm 
correspond to the two poling periods used in this letter, while 
the gray band indicates nonlinear-optical mirror behavior. 



In an inhomogeneous medium with some spatial peri- 
odicity, however, the phase matching condition is rewrit- 
ten to include the corresponding reciprocal lattice vector. 
In periodically poled QPM, for example, the sign of the 
nonlinear susceptibility is reversed in neighboring layers 
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with a period of A. Quasiphase matching then requires 

± k z (2u) = 2k z (u) + 2ra7r/A, (5) 

or, equivalent ly, A = mL co h, where m is any inte- 
ger. As mentioned earlier, poling in our system can be 
achieved by simple reorientation of appropriate sections 
of VLSSRs, analogous to domain engineering. Thus, tak- 
ing m = 1 in Fig. [3j the intersection of the dashed lines 
(representing poling periods of 12 and 24 cm) with the 
coherence length curves indicate the frequencies at which 
Eq. ([5| is expected to be satisfied, resulting in enhanced 
SHG output for the respective configuration. 

Furthermore, we see that something interesting hap- 
pens as the FF approaches the cutoff frequency. Near 
this frequency, k z (uj) vanishes. Thus, Eq. ^ can be 
satisfied for both the reflected and transmitted modes by 
the same poling period. This is indicated in Fig. [3] as the 
point where the reflected and transmitted modes cross, 
corresponding to where the FF passes through the cutoff 
frequency. Thus, for poling periods in the proximity of 
this coherence length, we expect to see QPM enhance- 
ment of both the reflected and transmitted SH waves. 

In order to inject and measure the FF and SH waves, 
respectively, the inner conductors (probes) of two coaxial 
cables are inserted into the waveguide on either side of the 
VLSRR sample. The probes are positioned 1.5 cm from 
the sample and covered with an alumina shell to reduce 
the impedance mismatch between the coaxial cable and 
the waveguide. Rectangular aluminum slabs with cross- 
section exactly equal to the waveguide are placed at a 
variable distance from the probes. By manually optimiz- 
ing the distance between these back walls and the probes, 
we are able to reduce reflections at the coaxial- waveguide 
interface to less than 0.75 dB at the SH frequencies of in- 
terest, preserving the directionality of the generated SH 
wave. The FF, however, does not need to be matched 
because the much larger material losses in the FF range 
prevent significant back-reflection inside the sample. In 
other words, the vast majority of the incident FF power 
that is able to enter the sample travels in the forward 
direction and is absorbed by the NLMM before reach- 
ing the second interface, preserving the directionality of 
the FF wave. Thus, by pumping through one probe at 
the FF and measuring the SH power at both probes, we 
are able to simultaneously investigate the reflected and 
transmitted SHG in the VLSRR loaded waveguide. For 
clarity, a diagram of the experimental apparatus is shown 
in Fig. [4] An Agilent PNA-X N5245A network analyzer 
is used as both the source and receiver. 

In order to validate our experimental results and in- 
terpretations, we employ nonlinear TMM calculations on 
a homogeneous slab with equivalent effective properties 
to the VLSRR loaded waveguide, sandwiched between 
semi-infinite slabs of vacuum, assuming the non-depleted 
pump approximation [I8j IT9] , The power incident on 
the sample at the FF is approximated by Pi = Pq — P r , 
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FIG. 4: (Color online) Diagram of the experimental setup 
employed to measure the SH spectrums, depicting a cross- 
section of the VLSRR loaded aluminum waveguide. 

where P r is the reflected power measured in the exper- 
iment and Po is the source's output power. The losses 
before and after the sample are assumed to be 0.75 dB at 
all frequencies. Although this system is clearly a rough 
approximation of the actual experimental setup, qualita- 
tive agreement is still expected, especially regarding the 
location and shape of the phase-matching enhancements. 

First, we investigate phase-matching of the reflected 
SHG in a negative-index band. All four sections are 
aligned to simulate a homogeneous nonlinear medium. 
The source is swept from 600 to 650 MHz with an output 
power of -5 dBm. The resulting transmitted and reflected 
SH spectrums are plotted in Fig. [5] The plot shows sig- 
nificant enhancement in the SHG of the reflected wave 
over the transmitted, and is supported qualitatively by 
the inset showing the TMM result. Furthermore, the 
peak frequency of 1255 MHz lies well within the range 
predicted by the coherence length calculation of Fig. [3] 
These results confirm the referenced theoretical studies, 
and in this configuration our NLMM can be considered 
the first realization of a nonlinear-optical mirror. More- 
over, for the network analyzer's maximum output, we 
were able to display nonlinear-optical mirror action with 
a conversion efficiency as high as 1.5%. 
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FIG. 5: (Color online) Comparison of the transmitted (solid 
blue) and reflected (dashed green) SHG powers when all sec- 
tions are aligned. The right inset depicts the experimental 
configuration, and the left inset shows the TMM calculation. 

In the next experiment, the previous configuration is 
retained in all ways except that sections 2 and 4 are phys- 
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ically rotated by 180 degrees, effectively producing a pe- 
riodically poled nonlinear crystal with A = 12 cm. The 
transmitted SH spectrum of this QPM setup is compared 
to that of the original in Fig. [6j The transmitted SH 
power is greatly increased due to QPM, with the peak 
frequency of 1253 MHz in excellent agreement with Fig. 
[3j Once again, qualitative agreement is found between 
the experimental data and the TMM calculations. 
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FIG. 6: (Color online) Comparison of the transmitted SHG 
power when the varactors are aligned (solid blue) and when 
they are periodically poled with A = 12 cm (dotted red). 
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FIG. 7: (Color online) Comparison of the transmitted (solid 
blue and dotted red) and reflected (dashed green and dotted 
black) SHG power when the varactors are aligned and when 
they are periodically poled with A = 24 cm, respectively. 



at microwave and terahertz frequencies. 
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In the final configuration, sections 1 and 2 are aligned 
and oriented opposite to sections 3 and 4, simulating a 
poling period of 24 cm. Since this poling period is close 
to the k z (uj) = crossing point in Fig. [3j we can ex- 
pect a strong enhancement of the reflected wave near a 
SH frequency of 1330 MHz, with a simultaneous (though 
weaker) enhancement of the transmitted wave. This is 
precisely what is shown in Fig. [7| with the SHG pow- 
ers for the un-poled case plotted for comparison. While 
the inset again shows good qualitative agreement be- 
tween the TMM and the experimental results, it should 
be noted that the peak emitted powers are off by an 
order of magnitude, whereas much better agreement is 
found in the previous two configurations. This is most 
likely due to the inherent sensitivity when approaching 
the near-zero frequency band, and it was confirmed in 
TMM calculations (not shown) that small changes in the 
linear properties produce large changes in the magnitude 
of the emitted radiation, but with similar qualitative be- 
havior. Moreover, it has been shown that spatial disper- 
sion, neglected here, can play an important role as the 
wavevector vanishes [22] . 

The exotic phase-matching effects shown here are made 
possible due to the unique linear properties accessible in 
NLMMs, in particular negative- and zero- index behavior. 
In addition, the use of varying periods to achieve peri- 
odic poling to switch between the various SHG schemes 
demonstrates the compatibility between NLMMs and 
QPM-a key result. We expect these results to be a nec- 
essary and vital step towards the realization of high effi- 
ciency nonlinear mirrors and other NLMM-based devices 
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